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Abstract 

Mechanics is developed over a differentiable manifold as space 
of possible positions. Time is considered to fill a one-dimensional 
Riemannian manifold, so having the metric as lapse. Then the sys- 
tem is quantized with covariant instead of partial derivatives in the 
Schrodinger operator. 

PACS-numbers: 0320 (Classical mechanics of discrete systems; general math- 
ematical aspects); 0365 (Quantum mechanics); 0240 (Differential geometry); 
9880 (Cosmology). 
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I. INTRODUCTION 



In General Relativity, the differential quotient between proper time and 
coordinate time is called lapse function. In the present article, this notion is 
used for an arbitrary classical mechanical system. Space is considered as n- 
dimensional Riemannian space V n and time is considered as 1-dimensional 
Riemannian space V\. Then the square of the lapse function turns out to 
be the metric of this V\. Possible applications and comparison with other 
approaches found in the literature will be shown in section VI below. 

Let us consider a mechanical system. The space of all possible positions 
shall be the n-dimensional differentiable manifold M n . It is endowed with 
local coordinates q l ,i = l,...n. Most of all mechanical systems have the 
property that M n is a subset of K m x.(S 1 ) n ~ m , so that the first m coordinates 
are Cartesian ones and the remaining are periodic ones (i.e., angles). Here, R 
denotes the space of reals, Z the space of integers, and the one- dimensional 
torus S* 1 can be defined as factor space S 1 = R/Z . But in general, M n cannot 
be covered by one single coordinate system. The time is denoted by t, and 
j- t will be denoted by a dot. So, <f is the velocity of a moving particle q l (t). 
Therefore, the velocity at time t is an element of the tangent space T x M n of 
M n at x = q l {t). The tangent bundle TM n is the union of all tangent spaces. 

Contrarily to the usual procedure we now introduce the lapse function 
N(t) which shall be an arbitrary positive function. (Here and below all 
functions shall have the necessary differentiability properties.) The proper 
time r is defined by 



It is uniquely determined up to an integration constant, i.e., without sped- 




(1) 
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fying the point where r = 0. The space of all possible times is a connected 
oriented one-dimensional Riemannian space V\ with coordinate x l = t and 
metric gn = N 2 (t). The orientation is chosen such that increasing time leads 
into the future. So, Eq. (1) represents the proper time r as proper length 
within this V\. 

Remark: The definition is chosen such that proper time does not depend 
on the velocity, so we do not cover relativistic effects. 

Each positive function N(t) defines a gauge, and results should not de- 
pend on it. In this manner, we define the following gauge-invariant quantity, 
the proper velocity v % 

v> = i i (2) 

We have to prove that v 1 does not depend on the special choice of N; this 
follows from Eqs. (1,2) via the equation 

i d <? 

The action I is the integral of a Lagrangian L 

I = J Ldt (3) 

and is supposed to be a coordinate-, gauge-, and T-invariant quantity. T- 
invariance means that / does not change if the orientation of V\ is reversed. 
The range of integration in Eq. (3) is a connected subset of Vi, i.e., any fixed 
time-interval; but we do not specify now which kind of interval is used. 
We restrict ourselves to first-order Lagrangians, i.e., L is a function 

L : TM n xV 1 — ► R (4) 

The next three steps are done by plausible arguments, not by proofs. 



First, the explicit t-dependence ( t e V\ ) of L, Eq. (4), is compatible 
with gauge-invariance of I only for the case that the t-dependence of L is 
via N(t) only, i.e., 

L = L(q\q\N) (5) 

Second, the coordinate- and gauge-invariance of / requires the following 
form of L 

L = G(q\v*)-N (6) 

where G is a certain scalar; this becomes plausible from Eqs. (1, 2, 3). 
Third, we assume that G can be developed into powers of v 1 

oo 

g = E4 fe U(?> n ---^ (7) 

k=0 

with certain tensors a^ k \ Here, and below, the Einstein sum convention is to 
be applied. Then it follows from T-invariance, that only even values k give 
a non-vanishing contribution to Eq. (7). 

The simplest non-trivial example for Eq. (7) is the case that only k — 
and k — 2 give contributions. To meet the usual notation we define 

V = h l3 = 2o4 2 V) (8) 

Inserting Eqs. (7, 8) into Eq. (6) we get 

L = 4//,,rV —V)-N (9) 

Without loss of generality, hij is assumed to be a symmetric tensor in M n . 
Here, the coordinate-, gauge-, and T-invariance of /, Eqs. (3, 9) is immedi- 
ately seen; so we also could have taken Eq. (9) as a definition of L. 
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To give the Lagrangian Eq. (9) the structure defined by Eq. (5) we insert 
Eq. (2) into Eq. (9) and get 

L = - 2 9ij <}' <}' - V-N (10) 



where we used the definition 



9a = ■ hij (11) 
Next, we introduce the momentum p^ by 

»-* < 12 » 

From Eq. (10) we get 

Pi = 9ij Q 1 (13) 

It holds: the momentum is gauge-invariant. This is proven by the fact that 
from Eqs. (2, 11, 13) one gets 

Pi = h i:j v j (14) 

From Eqs. (12, 13) we get 

d 2 L , , 

9ij (15) 



dcfdqi 

where gij depends on q % and N only. The analogous gauge-invariant equation 
to Eq. (15) reads 

w ^ h -> (16) 

and depends on q % only. 

Remark: One could use Eqs. (12, 15) also for the general case L, Eqs. 
(6, 7); but then would in general depend on the velocities, too. If g^ 
is interpreted as metric, then this would be the step from Riemannian to 



Finslerian geometry. A typical example of Finslerian geometry appears, if 
the term with k = 4 in Eq. (7) is allowed to appear. 
Let us introduce the Hamiltonian 



H = pa? - L (17) 

The canonical equations make sense only for the case that the velocities can 
be expressed as functions of the coordinates, momenta, and time. Looking 
at Eq. (13) one can see that this takes place if and only if g^ is a regular 
matrix. So, we assume this to be the case in the following and denote the 
inverse matrix to gij by g l K From Eq. (11) it follows that also hij is invertible. 
The inverse matrix to is denoted by h tj . It holds 

g ij = N .h i i (18) 

From Eq. (13) we get 

q l = 9 ij Pi (19) 
We insert Eqs. (10, 18, 19) into Eq. (17) and get 

h = \ g 13 PlPj + V-N (20) 

which can also be written as H = (| h 1 * Pi Pj + V) ■ N . The canonical 
equations are 

dpi 



and 



Eq. (21) is equivalent to Eq. (19), whereas Eq. (22) represents the equation 
of motion; in the next section we discuss it in more details. 
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II. THE EQUATION OF MOTION 

The acceleration is a 1 = q l . In general, the equation of motion expresses 
the acceleration as function of coordinates, velocity, and time. To get this 
structure, we insert Eqs. (13, 20) into Eq. (22). After some calculus we get 

a* = ^ € ~ V*- N - # q k Y) k (23) 

where V' 1 — g lj Vj and T denotes the Christoffel affinity (which is the same 
both for gij and hij). As usual, «,i » is an abbreviation for the partial 
derivative with respect to the coordinate << q l ». 

We can give three results immediately: First, for N and V being con- 
stant, the equation of motion is just the geodesic equation in the M n with 
Riemannian metric g^. Second, for iV and being constant, the equation 
of motion reads = a 1 + V' 1 and equals the classical equation of motion in 
the potential V. Third, using gauge-invariant quantities, we can write the 
equation of motion as 

0=— + r jk v^v k + h^V j (24) 

The first two terms of the r.h.s. represent the covariant derivative of the 
proper velocity with respect to proper time. 

In the next step we consider, independently of the Hamiltonian, under 
which condition the action I Eq. (3) has a stationary value. One should 
expect that the same equation of motion appears, but this is not fully trivial 
to show. 

The corresponding Euler-Lagrange equation to the action I reads 



With Eq. (12) we get 

A = | <») 
Comparing with Eq. (22) we have to show that 

™ = - S ± (27) 
dq l dq l V ; 

Looking at Eq. (17) one could get the impression that Eq. (27) can be 
fulfilled for a constant product piq 1 only, but this impression is wrong, because 
in the l.h.s., H is a function H(q\pi, N) but in the r.h.s., L is a function 
L(q\q\N). And so, with H Eq. (20) and L Eq. (10), the validity of Eq. 
(27) can be proven. 

III. THE LOWER-DIMENSIONAL CASES 

Let us consider the simplifications for the lower-dimensional cases. For 
n = 1, one knows that the Riemannian space V\ is flat, and so the Lagrangian 
Eq. (9) reduces to L = [^v 2 - V(x)] ■ N(t) with q 1 = x, v 1 = v and 
fin — rn — const. ^ 0. With N — 1 this is the usual point particle in a 
potential V. 

For n — 2, the Riemannian space V 2 = (M 2 , hij) need not to be flat, but 
it is always conformally flat. So one can always find local coordinates such 
that the Lagrangian Eq. (9) can be written as 

m Af 

L = l-v 2 + —w 2 - W(x, y)] ■ S(x, y) ■ N(t) (28) 

with q 2 = y,v 2 = w and /i 22 = M = const. ^ and W ■ S = V as additional 
relations. S ^ is the suitably chosen conformal factor. 

For n > 3, however, a V n need not to be conformally flat, and so, in 
general, the usual kinetic term with constant masses can be reached neither 
by a coordinate nor by a conformal transformation. 



IV. QUANTIZATION 

The usual quantization procedure is to substitute pk by %h-^ in the 
Hamiltonian to come from the function to the operator. If we make this 
in our approach, then gauge-invariance is automatically ensured, because 
both q k and pk are gauge-invariant quantities. (To prevent misunderstand- 
ings, we explicitly say: i is an index G {1, . . .n} if written in index position, 
and it is the imaginary unit otherwise.) But to ensure coordinate-invariance, 
the partial derivative is not sufficient. The most natural way to circumvent 
this difficulty is to use the covariant derivative with the same T as before. 
Then Vfc denotes the covariant derivative with respect to q k . 

The world function is denoted by ip, it is a function 

i> ■■ M n — > C (29) 

where C denotes the set of complex numbers. 

The energy of the system is E = H/N. It is a gauge-invariant scalar, 
and it is constant along classical trajectories: ^ = which follows from Eqs. 
(20, 21, 22). 

So we get the Schrodinger equation Hip = E ■ N ■ ip with ip — ip(q l ) and 

H = ~h 2 g ij V i V j + V-N (30) 

The zero energy Schrodinger equation simply reads 

% 2 Uip = 2V*p (31) 

where □ denotes the D'Alembertian with respect to the metric hij, i.e., □ = 
h lj ViVj, whereas the general Schrodinger equation can be obtained from this 
one by a suitable redefinition of V. 
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To circumvent the explicit calculation of the Christoffel affinities we apply 
the following formula 

□ = di Vh h ij dj (32) 

where h = |det/iy| ^ 0. 

Remark: One should observe that the form used here is surely the simplest 
possible way to get a coordinate-invariant Schrodinger equation; however, it 
is not the only possible one which goes over to the classical Schrodinger equa- 
tion (i.e., that one with partial derivatives) if becomes constant. Indeed, 
one could use the conformally invariant operator O c — D — £R instead of 

where R is the curvature scalar of the metric and £ = ^0^) ■ Only 
for n < 2 one has D c = □; for n = 2 because of £ = 0, and for n = 1 
because of R = 0. But even for n > 3 one can cover this variant by a suitable 
redefinition ofV. 

Let us shortly say what happens for the lower-dimensional cases. For 
n = 1, one simply uses coordinates such that hu = 1 and one gets the usual 
equation. For n — 2, however, it is a little more involved. We employ the fact 
that hij is conformally flat and so it can be written as h^ = y/h rjij where 
rjij is a matrix in diagonal form where all diagonal elements are G {+1,-1}. 
rfi is the inverse to 77^; and, by construction, they coincide. Then we insert 
Eq. (32) into Eq. (31) and get 

h 2 rf j didjip = 2 Vh V if; (33) 

The l.h.s. represents the flat -space D'Alembertian, and the factor y/h in the 
r.h.s. can be absorbed by a redefinition of V. 

For n > 3, however, it requires special circumstances to get the Schrodin- 
ger equation in the form of a flat-space D'Alembertian. 
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V. SOLUTIONS OF THE SCHRODINGER EQUA- 
TION 



From the full set of solutions of the Schrodinger equation (31) we are 
essentially interested in those solutions which correspond to the classical so- 
lutions of the system (21, 22). To this end we apply the WKB-approximation 
and insert the ansatz 

ij) = a- exp(iS/h) (34) 

into Eq. (31) and get 

h 2 Ua + ih(2a, k S' k + aDS) - aS, k S' k = 2aV (35) 

where S' k = h jk Sj. From Eq. (34) we have the situation that now two 
functions ( a, S ) represent one function ( ip ). So we are free to put an 
additional relation as calibration. It turns out that the following calibration 
is useful: we set for a moment h — 0, insert this into Eq. (35) and use the 
resulting equation 

S tk S' k + 2V = (36) 

as natural calibration. This is the usual classical limit. 

Before we proceed we must be sure that Eq. (36) possesses solutions. 
If the metric has indefinit signature, then this is trivial. Let be of 
definit signature; without loss of generality it shall be positively definit, for, 
otherwise, simply V has to change its sign. In regions where V < 0, Eq. 
(36) has solutions, but in regions with V > it does not have any solutions. 
One should remember here, that we have redefined V such that the whole 
system has zero energy. So, V > corresponds to a negative kinetic energy; 
the latter is impossible for a positively definit metric hij. We get as result: 
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the calibration Eq. (36) is possible if and only if classical motion takes place 
there. 

Now we insert Eq. (36) into Eq. (35) and get 

= h Oa + 2 i a k S' k + i a OS (37) 

To proceed, there exist different possibilities: first, one again neglects the 
term with h, second, one requires a to be a slowly varying amplitude such 
that Da is negligible in comparison with US, or, third, one thinks of a and 
S as real functions and so Eq. (37) splits into real and imaginary parts. It 
is not so essential which of these three arguments are applied, because all of 
them give rise to the equation 

= 2a k S' k + aDS (38) 

Eq. (38) can be solved as follows: let S(q l ) be a solution of Eq. (36) with 
S t k 7^ 0. There exists no time in the system, but we can introduce a time T 
by requiring that ^ = S' k d k . With b = In a 2 , Eq. (38) now reads 

% = -US (39, 

which can be integrated along the trajectories of T. In an afterwards- 
interpretation one can identify T with r, S ,k with v k and with p k ; this 
turns out to be compatible with the classical (= non-quantum) equations. 
But this alone does not suffice: from Eq. (39) one calculates the function 
a(g*) and inserts it together with S(q*) into Eq. (37). Then the WKB- 
approximation turns out to yield results close to the exact solution only for 
the case that indeed, \hOa\ is negligible in comparison to |aDS'|. So one can 
check in which region the semiclassical approach makes sense. 
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VI. CONCLUSION 

Classical mechanics, as is usually presented, e.g. in Refs. 1 and 2, uses 
essentially vector spaces as space of possible positions. Then one has the 
duality between coordinates and momenta (which we do loose here) and can 
build a symplectic manifold. Furthermore, one has usually a constant mass 
tensor (which means a constant matrix in our notation). 

Both points are generalized in the present paper. The present approach is 
inspired by work on Hamiltonian quantum cosmology, e.g. Ref. 3, where the 
space of possible positions is the set of all possible spatial geometries (called 
superspace). The set of all possible spatial geometries turns out to be neither 
a vector space nor is the matrix hij a constant one. Even, if one restricts to 
the minisuperspace which corresponds to homogeneous spatial geometries, 
one does not get a vector space. Example: The set of all homogeneous 

3- spaces of Bianchi-type IX [i.e., there exists a transitive subgroup of the 
isometry group isomorphic to 5*0(3)] which is a manifold with boundary, 
the interior is composed of points corresponding to spaces whose isometry 
group is 3-dimensional, and the boundary points are formed by spaces with 

4- dimensional isometry group (i.e., the axially symmetric Bianchi-type IX 
models), and the edge (the boundary of the boundary) consists of one line 
which itself corresponds to the isotropic 3-spheres with 6-dimensional isom- 
etry group. (Concerning details to this point see e.g. Ref. 4). 

If M n is such a manifold with boundary, then a trajectory is simply mir- 
rored at the boundary. 

Here we carefully distinguish between co- and contravariant tensor in- 
dices, and the Einstein sum convention is used in its strong version: summa- 
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tion over double indices takes place only for the case that one of them is in 
upper (= contravariant) and the other in lower (= covariant) position. It is 
a nice additional check of the formulae that the necessity to write the S-sign 
never appeared. 

The essential result of the present paper is to show that up to dimension 
two, the Schrodinger equation comes out with the flat-space D'Alembertian 
whereas for higher dimensions, it requires a special structure of the action to 
have this property. This has the following consequence for quantum cosmol- 
ogy: All models with one- or two-dimensional minisuperspace can be writ- 
ten with the flat-space D'Alembertian in the Schrodinger equation (which is 
called Wheeler de Witt equation here), whereas for higher-dimensional mini- 
superspace models, e.g. Ref. 5, this property requires a special structure of 
the underlying system. 

Reformulated for the classical (i.e., non-quantized) system one can state: 
A system with one or two degrees of freedom has always a kinetic energy 
which can be written as sum of terms of the type ± y v 2 with positive 
constant values m, whereas for three or higher dimension this need not to be 
the case. 

The kind of introducing the covariant derivative in Eq. (30) instead of the 
partial one is the mathematical background of the (today widely accepted) 
solution of the so-called factor-ordering problem, which filled many papers 
on quantum cosmology in the eighties, see Ref. 6 which is a bibliography of 
papers on the topic. 

We always wrote velocities with upper (contravariant) and momenta with 
lower (covariant) index; this is more than a purely notational arbitrariness, 
moreover, it is the only adequate form from the differential geometric point 
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of view. 

A geometric description of non-relativistic quantum mechanics has al- 
ready carried out by Kuchaf [7] in 1980. He uses a degenerate metric (i.e., a 
metric with vanishing determinant), so that he needs additional considera- 
tions to relate the co- and the contravariant components of it. He solves the 
factor-ordering problem by writing the Laplacian covariant with repect to 
this degenerate metric. Contrary to our approach (see also Ref. [5] for more 
details), he uses Dirac's constraint quantization. 

Section 7.2 of Ref. [8] develops classical mechanics in parametrized form. 
In this form, it becomes time-reparametrization invariant just as General 
Relativity is coordinate-invariant. Their approach takes velocities and mo- 
menta on the same footing (both are covariantly written vectors). 

The book [9] by Zeh reviews many aspects of the direction of time. In 
subsection 5.2.1 of that book, also the reparametrization invariance of time 
is mentioned, Zeh relates this property to Mach's principle (regarding time). 
Ref. [10] presents a geometrization of classical mechanics by use of a sym- 
plectic structure. Refs. [11] discuss the recovering of time and the deduction 
of the Wheeler de Witt equation in quantum cosmology. 

An application of the present approach (the present article is a revised 
version of the unpublished Potsdam-Report No. 93/10 from January 1993) 
can be found in section V A of Ref. [12], where it is used to deduce the 
Wheeler de Witt equation for the Starobinsky cosmological model. For fur- 
ther generalizations see e.g. Ref. [13]. 
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